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The free energy of a two-dimensional system at criticality has in general an universal part propor- 
tional the logarithm of the system size. This term was shown by Cardy and Peschel to be related to 
the curvature of the system, with smooth metrics and singular points contributing in distinct ways. 
In this paper we present a numerical study of the effect, for the Ising model on lattices with various 
topologies from the sphere to genus two surfaces. The nature of this term for specific lattice models 
is an open problem because the distinction between the two kinds of contributions involves an inter- 
change of the limit of singular curvature with the thermodynamic limit. For all the lattices studied 
we found precise agreement with the conformal field theory prediction for conical singularities. 



INTRODUCTION 

The free energy of a two-dimensional system, of char- 
acteristic size L, is expected from general arguments p| 
to have, at criticality and for fixed shape, a large L ex- 
pansion of the form 

F = .foL 2 + f b L + ClnL + D + o{l). (1) 

While the bulk free energy f and the boundary free 
energy /& depend on the system details the coefficients 
C and D are universal. The scale independent term D 
is known for strip geometries to be simply related to the 
theory central charge Q and for lattices without bound- 
aries to be a modular invariant quantity determined 
by the lattice shape. 

The coefficient of the logarithmic term C, from the 
conformal field theory point of view, was shown by Cardy 
and Peschel Q to be, for a system with a smooth metric 
and boundary, 

C = -\c X (2) 

where \ is the Euler characteristic of the system and c 
the central charge. This result can be understood as the 
integral, over the all system, of the trace anomaly which 
is proportional to the local curvature. However in the 
presence of scale invariant geometric shapes, as bound- 
ary wedges or a bulk conical singularities, the expression 
above fails to hold. For instance, for a conical singularity 
with deficit angle e = 2ir — 9, where 9 is the angle spanned 
by the cone, the contribution to the free energy is Q] 

C e = ^(l-{2*/9f) (3) 

and not a term proportional to —e as one would expect 
from a delta function singularity on the curvature. 
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This result is surprising because it signals that, con- 
cerning the free energy, a conical singularity cannot be 
considered as a limiting case of a smooth metric. Ap- 
parently the thermodynamic limit, with the ultraviolet 
cutoff going to zero, does not commute with the limit of 
taking the curvature at a point to infinity. It is then natu- 
ral to ask, following reference Q , which of the expressions 
© or ©, if any, describes the logarithmic correction C 
for a model defined on a regular lattice. 

In this paper we consider this problem for the Ising 
model defined on square lattices, without boundaries, 
on various topologies from the sphere to genus two sur- 
faces. The Ising model on lattices with the topology of 
the sphere has been studied previously but the coef- 
ficient C was obtained only for a restricted class of hexag- 
onal lattices 0. 

We consider three lattices with the topology of the 
sphere; the cube lattice, its dual and the L-shaped lattice 
shown in Fig. ^ We consider also the genus two lattice 
shown in Fig. [21 that can be visualized as a torus with 
an additional handle All these lattices are locally 
fiat, equivalent to the regular square lattice, with the 
exception of a few corners that can be seen as conical 
singularities. For the cube and the L-shaped lattice these 
are the vertices where three or five squared faces meet. In 
the dual cube lattice the singularities are located on the 
triangular faces, that share edges with only three squared 
faces by opposition to four. For the genus two lattice of 
Fig. [21 the singularities correspond to the two octagonal 
faces, shown in dashed line, that have eight neighboring 
squared faces. 

Closed form expressions of the Ising model free en- 
ergy are known for toroidal ||, cylindrical and some 
non-orientable [ToL Hl| geometries, where translational 
invariance permits the use of Fourier transforms. For 
the lattices considered in this paper, such an analytical 
treatment is not possible. The free energy will be evalu- 
ate numerically, using the Kasteleyn method, for a series 
of lattices with fixed shape and increasing size. 



The logarithmic correction C that we obtain for the 
various lattices considered, corresponds very accurately 
to the sum over the lattice singularities of the conical 
singularity contributions Cg as given by Eq. l|5Jl. 



EVALUATION OF THE FREE ENERGY 
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The Ising model partition function, on a two dimen- 
sional lattice, can be expressed using the Kasteleyn 
method 0, Q, E| in terms of the Pfaffians of certain 
antisymmetric adjacency matrices. The formalism has 
a combinatorial interpretation in terms of the high tem- 
perature polygon expansion of partition function. It can 
alternatively by seen as graphical prescription to express 
the Ising model partition function as a Grassman inte- 
gral. The resulting Pfaffians, for an arbitrary lattice, 
cannot in general be evaluated in closed form but pro- 
vide a good numerical tool to calculate the free energy 
for lattices with reasonably large size. In this section we 
will review general aspects of the Kasteleyn formalism 
needed to evaluate the partition function for the lattices 
of Fig. ^ and [21 For additional details andproofs of the 
method we direct the reader to references 0, EI IHI ■ 

Consider a lattice, or more properly a graph, G drawn 
without superposition of edges on a surface of genus g, 
where the genus is the number of handles in the sur- 
face related with the Euler characteristic by % = 2 — 2g. 
The Ising model on such a lattice is defined by assigning 
a coupling constant to each edge connecting two lattice 
vertices and placing an Ising spin at each vertex. 

Following the Kasteleyn method we start by decorat- 
ing the lattice vertices with a decoration graph, meaning 
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FIG. 2: The higher genus lattices: a) the usual Mi rows by 
Mi columns toroidal lattice and b) a genus two lattice char- 
acterized by five integer sizes (Mi : i = 1, . . . , 5). The identifi- 
cations of edges at the boundaries is given by the letters. The 
genus two lattice can be seen as a torus with an additional 
handle on the bulk. 
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FIG. 1: The lattices with the topology of the sphere: the 
cubic shaped lattice, its dual and the L-shaped lattice. Each 
lattice is characterized by a set of integer dimensions Mi. 



that each vertex on the lattice G is replaced by a cluster 
of vertices. In the Grassman formulation, each vertex 
of the decoration corresponds to a Grassman variable 
associated with the original lattice site. The resulting 
decorated lattice will be denoted by Gd- Our choice of 
decoration for vertices with coordination number three, 
four and five is shown in Fig. [3] The order by which the 
exterior edges connect to the decoration polygon is not 
important and the decoration graphs can be rotated. 

The edges of the decorated lattice Gd are assigned 
weights: 1 to the internal edges of the decorations and 
a weight w = tanhi-T to the edges inherited from the 
lattice G, with K being the Ising coupling constant on 
that edge in units of KbT. The edges of Gd need also to 
be oriented, assigned a direction that can be represented 
graphically by an arrow. The internal edges of the deco- 
rations have already an orientation given in Fig. [3J The 
remaining edges, inherited from the original lattice G, 
are oriented according to the Kasteleyn rule: in such a 
way that all lattice faces have an odd number of clock- 
wise oriented edges. This condition can be understood 
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Schematically 
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z if i — > j 
-z if j i 
otherwise 



(4) 





FIG. 3: The vertex decorations used for vertices with three, 
four and five incoming edges. 



as a consistency condition for the Grassman integral de- 
scription, corresponding to the concept of spin structure 

in a spinor formalism. 

An example of a decorated lattice orientation, that sat- 
isfies the Kasteleyn condition, is shown in Fig. for the 
cube lattice with dimensions (Mi, M2, M3) — (4,5,4). 
For a given lattice Gd there are many different edge ori- 
entations that verify this condition. In fact given a Kaste- 
leyn orientation we can generate a distinct Kasteleyn ori- 
entation by simply reversing the orientations (or arrows) 
on all the edges meeting at a given vertex. Two lattice 
orientations related by such an arrow reversal, or a series 
of arrow reversals, are said to be equivalent. 

It can be showed 0, ITsL [fif that for a genus g lat- 
tice there are precisely A 9 un-equivalcnt Kasteleyn orien- 
tations, matching the number of spin structures on the 
corresponding continuum description. This number is de- 
termined by the dimensions of the first homology group 
of the embedding surface. If aj, bi with i = 1, . . . ,g are 
loops in the surface forming a canonical basis of the first 
homology group then the A 9 un-equivalent orientation 
can be generated from an initial Kasteleyn orientation 
by reversing the orientation of the edges crossed by the 
a choice of such loops. 

Given a Kasteleyn orientation of the decorated graph 
Gd and a labeling of its vertices with an integer from 1 to 
A4, we define the corresponding adjacency matrix A as 
a M. x M. matrix with entries Aij that vanish if vertices 

1 and j are not connected by an edge, and take a value 
±2 if vertices i and j are connected by an edge of weight 
z, the sign being determined by the edge orientation. 



The Ising model partition function on the lattice G can 
then be expressed as 



Z Isins (K) = 2 Nv (cosh K) 



I 



4" 

- 5>,Pf MK) (5) 

i=l 



were Ny and Ne are the number of vertices and edges in 
G and we assume that all edges have the same coupling 
constant K . The sum runs over representatives of the A 9 
un-equivalent Kasteleyn orientations. The on take values 
±1 and are completely determined by the arrow parity of 
the non-trivial topology loops ai,bi along lattice edges. 
Modulo these signs the Pfaffians above do not depend on 
which orientation is chosen from each equivalence class 
or on the choice of the labeling of vertices. 

For translational invariant lattices, such as the toroidal 
square lattice, the Pfaffians of these adjacency matrices 
can be evaluated in a closed form. For the lattices shown 
in Figs. ^and[21such an analytic treatment is not possible 
and we are forced to resort to numerical evaluations of 
the Pf Ai , the Pfaffian of an antisymmetric matrix being 
the square root of its determinant. 
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FIG. 4: An orientation of the decorated cube lattice that 
verifies the Kasteleyn rule. The edges along the boundary are 
identified to form a cube with dimension 4x5x4. At each 
lattice site there is a four point decoration as shown in the 
inset. 
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LATTICES WITH THE TOPOLOGY OF THE 
SPHERE 

For the genus zero lattices shown in Fig. ^ the par- 
tition function is expressed, according to Eq. (J5J, by 
the Pfaffian of a single adjacency matrix. This corre- 
sponds to single Pfaffian used by Kasteleyn to evalu- 
ate dimer coverings on the rectangular lattice with edges, 
notice that a spherical lattice can always be deformed and 
drawn in the plane without superposition of edges. 

We evaluated the free energy for sequences of lattices 
with fixed shape and increasing size, by taking the integer 
dimensions characterizing the lattice Mj to be of the form 
Mi = rrii L with fixed rrii and increasing L. The Ising 
coupling constant was fixed for all edges at the square 
lattice isotropic critical value 



sinh2i4T r = 1. 



(0) 



Table IIVI shows typical examples of the free energies 
obtained, in units of KbT. The dual of the cube lattice 
was considered as a test to the stability of the numerical 
algorithms for large lattice sizes. The partition function 
on a given planar lattice Z and on the corresponding dual 
lattice Z* are related, at the critical point JBJ, by 

\ogZ{K c ) = log2(AV _ _E _ i) + log Z*(K C ) (7) 

with (Ny — — 1) = 1 for the cube lattice. This con- 
dition is a non-trivial check on the precision of the nu- 
merical evaluation of the Pfaffians in Eq. (0. From the 
values given in table IIVI we see that the condition @ 
is verified with a precision better than 1 in 10 9 and we 
can trust that the evaluation of the determinants of the 
adjacency matrices is not impaired by the large lattice 
sizes. 

The logarithmic correction C can be found by sub- 
tracting the leading volume term fo to the free energy, 



lattice: 


cube 


cube-corner 


L-shaped 


shape: 


(111) (121) 


(111) 


(11111) (12111) 


fits 


-0.19468 -0.19043 


-0.18675 


-0.20531 -0.19827 




-0.19465 -0.19100 


-0.18938 


-0.20566 -0.20150 




-0.19463 -0.19142 


-0.19093 


-0.20573 -0.20275 




-0.19461 -0.19175 


-0.19190 


-0.20572 -0.20338 




-0.19459 -0.19201 


-0.19254 


-0.20571 -0.20376 




-0.19458 -0.19222 


-0.19298 


-0.20569 -0.20402 




-0.19456 -0.19240 


-0.19329 


-0.20567 


expected: 


-0.19444 


-0.20556 



TABLE I: Successive fits to the logarithm correction C for 
lattices with the topology of the sphere and various shapes 
(Mi, . . . , Mat). Each value corresponds to the slope obtained 
by a linear regression using four consecutive free energy val- 
ues. There is a smooth convergence to the predicted correc- 
tion with increasing lattice size. 
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FIG. 5: The logarithmic correction for lattices with the topol- 
ogy of the sphere: the value of the residual free energy is given 
in function of the logarithm of the number of lattice vertices 
Ny. The uncertainties on the values are smaller than the 
symbols size. 



that for a square lattice is [Lj 



/o = 2G/7r+-ln2 



(8) 



where G is the Catalan constant. 

In Fig. |S]the residual free energy F—f Ny is plotted as 
function of the logarithm of the lattice size, log Ny /2, for 
various lattice shapes. We find a simple linear behavior 
and the coefficient C is the slope of the curve in the limit 
of large lattice size. The factor of 1/2 follows from QJ, we 
want the the logarithmic correction on the length scale 
y/Ny an d n °t on the lattice area Ny. 

We evaluate a series of values converging to C by do- 
ing linear regressions on sets of four consecutive lattices 
with increasing sizes. The values obtained for the three 
lattices, with various shapes, are shown in table [H We 
find a clear convergence to a value of C, that is the sum 
over the lattice conical singularities of the Cg correction 
inEq. (0. The rate of convergence depends on the shape 
or aspect-ratios of the lattice, being fastest for the more 
symmetric lattices. 

The conical singularities are of two kinds: vertex sin- 
gularities for the cube and L-shaped lattices and face 
singularities for the dual cube lattice. This distinction 
is artificial since the two types are related by a duality 
transformation. We can therefore consider only vertex 
singularities, with three and five square faces meeting at 
a single vertex, corresponding respectively to 6 angles 
of 37r/2 and 57r/2. The corresponding correction Cg are 
listed on Table HT1 The expected logarithmic correction 
C for a given lattice, see tabled is obtained by summing 
the Cg for all singularities occurring in that lattice. 

The nature of the conical singularities demands fur- 
ther clarification, since it is not immediately clear how 
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FIG. 6: The cube lattice with a corner removed. The new 
corner has six vertices where three or five squared faces meet 
but the conical singularity remains the same. Regardless of 
the detailed corner structure, both lattices at large scale fold 
into a 8 — 3ir/2 conical shape around each corner. 



to related the combinatorics of the lattice structure with 
continuum concepts as conical singularities in the metric 
of a surface. For the purposes considered in this paper, 
namely finite size corrections to the free energy, the im- 
portant aspect to obtain a conical singularity behavior 
seems to be the large scale structure of the lattice and 
not the detailed small scale structure. 

Consider for instance the cube lattice with a corner 
removed, as shown in Fig. In the cube lattice three 
squared faces meet at each corner. In the lattice of Fig. 
one of the corners has three vertices where three squared 
faces meet and three vertices where five squared faces 
meet. Regardless of these local differences, at large scale 
the corner for both lattices can be seen as a sheet of flat 
square lattice folded into a 9 — 3tt/2 conical shape. It is 
this fact that determines the conical singularity behavior 
and not the detailed lattice combinatorics. 

That this is the case can be seen from Fig. [5]and table 
H] where it is shown that the cube lattice with the corner 
removed has the same logarithmic correction C as the 
simple cube lattice. The different rate of convergence 
can be seen as due to the lattice defects introduced by 
the changes in the corner. 



HIGHER GENUS LATTICES: THE TORUS AND 
GENUS TWO 

We now consider the genus two lattice shown in Fig. |3 
the toroidal lattice, whose logarithmic correction is zero, 



# squares 


e 6 


Ce 


3 


n/2 3tt/2 


-0.024306 


5 


-tt/2 5tt/2 


0.01875 


8 


— 2TY 47T 


0.0625 



TABLE II: The deficit angle e, the spanning angle and 
the corresponding Cg correction for conical singularities with 
three, five and eight squared faces meeting at a vertex or face. 



lattice: 


torus (If) 


genus 2 (11111) 


Pfafhan: 


(+ -) (+ +) 


(-- + -)( ) 


fits 


-0.00020 0.00080 
-0.00017 0.00070 
-0.00016 0.00062 
-0.00014 0.00056 
-0.00012 0.00050 
-0.00011 0.00045 
-0.00010 0.00041 


0.12858 0.12942 
0.12819 0.12895 
0.12786 0.12855 
0.12759 0.12822 
0.12735 0.12793 
0.12715 0.12768 
0.12698 0.12747 


expected: 





0.125 



TABLE III: Successive fits to the logarithm correction C for 
the toroidal and genus two lattices with shapes (Mi, . . . , Mn), 
using different Kasteleyn orientations. The Kasteleyn orien- 
tations are denoted by ± signs labeling the parity along the 
lattice cycles with non-trivial topology. 



will also be studied as a non-trivial check of the methods 
used. 

For a genus g lattice the partition function is expressed 
in terms of the Pfaffians of 4 9 adjacency matrices, four 
for the torus and sixteen for the genus two lattice. Some 
of these Pfafhan vanish at criticality in the thermody- 
namic limit, while the remaining ones will give the same 
contribution to the logarithmic correction. This follows 
from the fact 0, that the various Pfaffians converge 
rapidly to a common bulk term times a topology and 
shape determined factor, 



Pf Ai(N v ) = ei Z(Nv) 



(9) 



where the &i are constants for large enough lattice size 
Ny. The free energy F follows from JSJ to be 



—F = N v In 2 + N e In cosh K 




InZ. 



(10) 



The logarithmic correction is entirely due to the Z 
term and can be equivalently evaluated from any non- 
vanishing Pfaffian. For this purpose it is convenient to 
introduce the auxiliary quantities 



F% (Ny) = F(N V ) - In 9* + In ( i £ a k Q,}j 



(11) 



where i is a label of the 4 9 un-equivalent Kasteleyn orien- 
tations. In the thermodynamic limit the Fi are distinct 
from the free energy F by a constant only. 

For both lattices we evaluate the determinants of the 
adjacency matrices by keeping the lattice aspect ratios rrii 
fixed and increasing the lattice size according to a scale L, 
with the lattice dimensions being Mi — rriiL. Examples 
of the evaluated Fi, for different Kasteleyn orientations, 
are shown in table HVl for various lattice sizes. 

Subtracting the leading volume contribution to the free 
energy, the residual free energy F — faNy is found to 
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FIG. 7: The logarithmic correction for the torus and the genus 
two lattice: the value of the residual free energy is plotted in 
function of the logarithm of the number of lattice vertices Nv ■ 
For the torus and the genus two lattice, respectively two and 
three un-equivalent adjacency matrices are considered. 



have a linear behavior on the logarithm of the lattice size 
logiVy/2, see Fig. [7| The coefficient C is the slope of 
these curves in the thermodynamic limit. As in the pre- 
vious section we generate a sequence of values converging 
to C by doing linear regressions on sets of four lattices 
with consecutive lattice sizes. The results obtained are 
shown in table Hill 

As expected the torus behaves as a flat lattice, while 
the genus two lattice reproduces the correction © ex- 
pected for two conical singularities with deficit angle 
e = —2ir. These singularities correspond to the two oc- 
tagonal lattice faces, shown in dashed line in Fig. [21 
that in a dual description are vertices where eight squares 
meet. 

The independence of the results on the choice of the 
Kasteleyn orientation is observed in Fig. [7] curves for 
different Fi show equal slopes. The vertical displacement 
is due to different log 9^ among the Pfaffians of the dif- 
ferent adjacency matrices. 

For the toroidal lattice the coefficient C vanishes at a 
rate consistent with the rates of convergence observed for 
C in the cube and L-shaped lattices. The convergence is 
somewhat slower for the genus two lattice, for which the 
rate is similar to the ones observed on the less symmetric 
genus zero lattices. For all cases the validity of Eq. © 
is established without doubt. 



CONCLUSIONS 

In this paper we have provided an answer to the ques- 
tion, posed by Cardy and Peschel in 0] , to which of the 
two forms © or for the logarithmic correction to the 
free energy is observed in specific lattice models. We 





— F for cube lattice 




— Fi for genus 2 lattice 


L 


(1,1,1) 


dual (1,1,1) 


L 


(--+-) 


( ) 


11 


559.558730 


558.865583 


14 


729.106047 


729.125841 


12 


676.718852 


676.025704 


15 


836.941813 


836.961544 


13 


805.033712 


804.340565 


16 


952.215747 


952.235424 


14 


944.503569 


943.810422 


17 


1074.92777 


1074.94740 


15 


1095.12862 


1094.43547 


18 


1205.07782 


1205.09741 


16 


1256.90902 


1256.21588 


19 


1342.66584 


1342.68540 


17 


1429.84491 


1429.15176 


20 


1487.69179 


1487.71132 


18 


1613.93638 


1613.24323 


21 


1640.15563 


1640.17514 


19 


1809.18352 


1808.49038 


22 


1800.05734 


1800.07682 


20 


2015.58641 


2014.89326 


23 


1967.39687 


1967.41633 



TABLE IV: Free energy F for the cube lattice and its dual 
with aspect ratios (Mi,M2,Ms) — (1,1,1) and the Fi free 
energies for the genus two lattice with aspect ratio(Affc : k = 
1, . . . , 5) = (1, 1, 1, 1, 1), for two distinct Kasteleyn orienta- 
tions. 



found that for regular lattices with conical singularities 
the logarithmic correction C is given by Eq. J3J). This 
was observed for a number of lattice boundary conditions, 
with various lattice topologies and shapes, provided that 
the lattice is locally a flat regular lattice with the excep- 
tion of a few vertices where at large scale it can be seen 
as folded into a conical shape. 

An interesting question that remains to be answered 
is if there are lattice realizations of Eq. J5J), the smooth 
metric contribution. 
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